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Abstract. The version 6 cloud products of the Atmospheric 
Infrared Sounder (AIRS) and Advanced Microwave Sound- 
ing Unit (AMSU) instrument suite are described. The cloud 
top temperature, pressure, and height and effective cloud 
fraction are now reported at the AIRS field-of-view (FOV) 
resolution. Significant improvements in cloud height assign- 
ment over version 5 are shown with FOV-scale comparisons 
to cloud vertical structure observed by the CloudSat 94 GHz 
radar and the Cloud-Aerosol LIdar with Orthogonal Polar- 
ization (CALIOP). Cloud thermodynamic phase (ice, liquid, 
and unknown phase), ice cloud effective diameter ( D c ). and 
ice cloud optical thickness (t) are derived using an optimal 
estimation methodology for AIRS FOVs, and global distri- 
butions for 2007 are presented. The largest values of r are 
found in the storm tracks and near convection in the tropics, 
while D e is largest on the equatorial side of the midlatitude 
storm tracks in both hemispheres, and lowest in tropical thin 
cirrus and the winter polar atmosphere. Over the Maritime 
Continent the diurnal variability of r is significantly larger 
than for the total cloud fraction, ice cloud frequency, and D e , 
and is anchored to the island archipelago morphology. Impor- 
tant differences are described between northern and southern 
hemispheric midlatitude cyclones using storm center com- 
posites. The infrared-based cloud retrievals of AIRS provide 
unique, decadal-scale and global observations of clouds over 


portions of the diurnal and annual cycles, and capture vari- 
ability within the mesoscale and synoptic scales at all lati- 
tudes. 


1 Introduction 

Clouds remain the largest source of uncertainty in future cli- 
mate projections (IPCC AR4; see Appendix A for list of 
acronyms). Several global and multidecadal observational 
data sets of cloud amount, cloud top height, optical thickness 

(r), effective radius (r e ), and cloud type are readily avail- 
able for addressing this source of uncertainty. These include 

(but are not limited to) the International Satellite Cloud Cli- 
matology Project (ISCCP; Rossow and Schiffer, 1999), the 
High-resolution Infrared Radiation Sounder (HIRS; Wylie 

et al., 2005), and the Advanced Very High Resolution Ra- 

diometer (AVHRR; Heidinger and Pavolonis, 2009; Foster 

and Heidinger, 2013). Over the last decade, an advanced set 

of complementary observations of cloud top properties and 
cloud vertical structure has been obtained with NASA’s A- 
train constellation of satellites. In particular, the Moderate 
Resolution Imaging Spectroradiometer (MOD1S; Platnick et 
al., 2003) has provided a wide array of 1 and 5 km resolu- 
tion cloud products from both the Terra and Aqua platforms 
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since 1999 and 2002, respectively. Furthermore, a much bet- 
ter understanding of the global vertical cloud structure has 
been obtained from the CloudSat 94 GHz radar (Stephens et 
al., 2008) and the Cloud-Aerosol LIdar with Orthogonal Po- 
larization (CALIOP; Winker et al., 2010). 

Ongoing satellite data set comparisons have revealed that 
discrepancies among an assortment of publicly available 
satellite cloud data sets can be explained by differences in 
instruments, algorithms, and sampling (e.g., Stubenrauch et 
al., 2013). These multidecadal data sets (and their associ- 
ated instrument simulators) have been invaluable for process- 
based evaluations of climate models (Lau and Crane, 1995; 
Klein and Jakob, 1999; Zhang et al., 2005; Williams and Tse- 
lioudis, 2007; Kay et al., 2012; Pincus et al., 2012). How- 
ever, their utility for long-term, global-scale cloud trends has 
been uncertain and difficult to determine (Wylie et al., 2005; 
Evan et al., 2007; Norris and Slingo, 2009). The sign and 
magnitude of a particular trend may strongly depend on the 
cloud type, geographical region, and geophysical parameter 
of interest (Dim et al., 2011; Bender et al., 2012), and satel- 
lite sampling characteristics may complicate the assessment 
of the diurnal cycle (e.g., Foster and Heidinger, 2013). The 
reasons for discrepancies among satellite estimates of cloud 
presence, amount, cloud top temperature, and effective emis- 
sivity are better understood (e.g., Rossow et al., 1985; Nasiri 
et al., 2011) than the differences among various estimates 
of r, ice cloud effective diameter ( D e ) (Stubenrauch et al., 
2013), and cloud thermodynamic phase (Chylek et al., 2006; 
Nasiri and Kahn, 2008; Hu et al., 2010; Riedi et al., 2010). 

Tropical trade cumulus (e.g., Medeiros et al., 2008) and 
subtropical stratocumulus (e.g., Bony and Dufresne, 2005) 
are key players in controlling climate sensitivity to radiative 
forcing. Recent studies have also shown the roles that cloud 
top height, cloud thermodynamic phase, and ice cloud mi- 
crophysics play in determining climate sensitivity. Zelinka 
et al. (2012) showed that the positive longwave cloud feed- 
back due to rising cloud height in the tropics and extratropics 
outweighs the negative cloud feedback from the reduction in 
high-cloud amount. Clement et al. (2009) described obser- 
vational evidence of a reduction in low-cloud amount with 
increasing SST that leads to a positive shortwave feedback. 
Other studies find evidence for a negative cloud feedback in 
the middle and high latitudes, and possibly from an increase 
in optical depth and/or a transition from ice to liquid phase 
rather than an increase in cloud amount (Gordon and Nor- 
ris, 2010; Zelinka et al., 2012). Trenberth and Fasullo (2010) 
correlated the realism of current-day Coupled Model Inter- 
comparison Project Phase 3 (CMIP3) simulations of South- 
ern Ocean subtropical cloudiness to global estimates of cli- 
mate sensitivity. This is a region with highly uncertain cloud 
characterization, including the spatial and temporal morphol- 
ogy of cloud thermodynamic phase. 

Recent climate model evaluations highlight an emerging 
need for additional observational constraints of ice cloud mi- 
crophysics (Hendricks et al., 2011; Gettelman et al., 2010, 


2012) and thermodynamic phase (Tsushima et al., 2006; 
Cheng et al., 2012), including the complex characteristics 
of mixed-phase clouds (Storelvmo et al., 2008; Klein et 
al., 2009). Many contemporary climate models contain ex- 
plicit representations of ice nucleation, ice supersaturation, 
and multiple types of cloud and precipitating hydrometeors. 
An example is the NCAR Community Atmospheric Model 
(CAM) version 5 (CAM5), which provides more realistic ice 
water content, cloud frequency, and mixed-phase cloud dis- 
tributions along with new physical representations of cloud 
processes. Climate models are sensitive to the formulation of 
their ice physics (e.g., Barahona et al., 2010; Gettelman et 
al., 2010; and references therein). 

Fasullo and Trenberth (2012) demonstrated that the CMIP 
Phase 5 (CMIP5) models that contain the lowest subtropical 
relative humidity (RH) minima are not only the most real- 
istic when compared to AIRS RH but also have the high- 
est climate sensitivities. The connections between climate 
sensitivity, temperature, and humidity, including their scale- 
dependence and temporal variability, cannot be fully under- 
stood independent of cloud morphology (Wood and Field, 
2011; Quaas, 2012; Fasullo and Trenberth, 2012). The phys- 
ical properties of clouds and the larger thermodynamic fields 
that they are embedded within must be linked at the instanta- 
neous temporal and FOV scales of existing satellite data sets. 

Numerous studies have demonstrated that the hyperspec- 
tral infrared observations offer greater potential to quan- 
tify ice cloud microphysics and cloud thermodynamic phase 
compared to broadband observations (e.g., Wei et al., 2004; 
Nasiri and Kalin, 2008; Guignard et al., 2012; Smith et al., 
2012; Wang et al., 2013). However, hyperspectral 1R cloud 
retrievals are difficult to calculate because of an excessive 
computational expense, the complexity and variability of 
cloud geometry, uncertainties in the underlying surface and 
atmospheric state, the necessity of using a large set of chan- 
nels (2378 between 3.7 and 15.4 pm in the case of AIRS), 
and the ongoing difficulty of developing an automated, rig- 
orous, and fast retrieval that is applicable to a 10+ yr data 
record with approximately 2.9 million AIRS FOVs per day. 
In the work described here, we adopt the optimal estima- 
tion retrieval approach of the Tropospheric Emission Spec- 
trometer (TES; Bowman et al., 2006; Kulawik et al., 2006; 
Eldering et al., 2008) as an a posteriori retrieval algorithm 
for the ice cloud parameters. This algorithm is computation- 
ally rapid and theoretically rigorous, and is further increased 
in computational efficiency by holding constant all surface 
and atmospheric parameters except for r, D e , and Tq in ice 
clouds. 

This paper describes improvements to existing cloud 
parameters, and an approach to retrieve new key cloud 
parameters from the Atmospheric Infrared Sounder 
(A1RS)/Advanced Microwave Sounding Unit (AMSU; 
Aumann et al., 2003; Chahine et al., 2006) that ultimately 
can be used in synergy with the retrieved temperature and 
specific humidity profiles in atmospheric investigations. The 
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Cross-track Infrared Sounder (CrIS) on the Suomi National 
Polar-orbiting Partnership (Suomi NPP), to be launched on 
the Joint Polar Satellite System (JPSS), and the Infrared 
Atmospheric Sounding Interferometer (IASI; Hilton et al., 
2012) on EUMETSAT’s MetOp satellite will continue the 
AIRS legacy into the future. 

The paper is organized as follows. In Sect. 2, we de- 
scribe the AIRS level 2 (L2) standard cloud products, im- 
provements made in the version 6 release, and demon- 
strate improved AIRS cloud top height retrievals when com- 
pared against CALIOP and CloudSat cloud top observations 
matched at the FOV scale. In Sect. 3, the AIRS L2 support 
cloud products are described including cloud thermodynamic 
phase and ice cloud r, 7) e , and Tq. The method of the ice 
cloud property retrieval, the logic behind and application of 
the quality control (QC), and an initial climatology for 2007 
are presented and described. In Sect. 4, the AIRS cloud prod- 
ucts are used in two separate “stress tests”: (1) diurnal vari- 
ability of ice cloud properties in the Maritime Continent, a 
region with a large convective diurnal cycle, and (2) compos- 
ites of summertime Southern Hemisphere (SH) and winter- 
time Northern Hemisphere (NH) midlatitude cyclone cloud 
properties. We discuss and summarize the major findings in 
Sect. 5 and describe a future outlook in Sect. 6. 


2 The AIRS instrument and cloud top 

pressure/temperature and effective cloud fraction 
retrieval 

The AIRS instrument suite was launched on the Aqua satel- 
lite in May 2002, and is in a sun-synchronous orbit with an 
equatorial crossing time of 01:30 and 13:30 local time (LT) in 
the descending and ascending node, respectively. The spatial 
resolution of AIRS at nadir view is 13.5 km and there are ap- 
proximately 2.9 million FOVs observed in a single day. The 
temperature and water vapor profiles are obtained in 3 x 3 ar- 
rays of AIRS FOVs that are coregistered to AMSU in a field 
of regard (FOR), and up to 324 000 retrievals per day are 
obtained by a cloud-clearing methodology (Susskind et al., 
2003, 2006; see Appendix B). Some of the major changes 
between the operational version 5 (v5) and version 6 (v6) re- 
trieval algorithms are summarized in Olsen et al. (2013). The 
cloud top height (Zc), pressure (Pc), temperature (7c), and 
effective cloud fraction (EOF) fields are retrieved after com- 
pletion of the cloud-clearing steps by comparing calculated 
and observed AIRS radiances in a set of channels sensitive 
to cloud amount and height (Kahn et al., 2007). The EOF 
retrieval assumes the cloud emissivity is spectrally flat. The 
EOF is the product of the cloud fraction and the cloud emis- 
sivity, i.e., the value of cloud fraction if the emissivity were 
always 1.0. For a cloud that is transmissive, the equivalent 
opaque fraction is reported, as this algorithm cannot distin- 
guish radiation through clouds from radiation around clouds. 


For simplicity, and the availability of correlative data sets, the 
cloud top 7c and Zc will be presented henceforth. 

2.1 What is new in version 6 

The most significant change between the v5 and v6 cloud 
retrieval algorithm is that Tq is retrieved within every nom- 
inal 13.5 km diameter AIRS FOV instead of within the en- 
tire 45 km diameter A1RS/AMSU FOR. Version 5 uses a sin- 
gle retrieval with 20 parameters: two layers of Tq for the 
AIRS/ AMSU FOR, and two layers of nine EOF values for 
each of the nine AIRS FOVs. Version 6 uses nine separate re- 
trievals in the AMSU FOR and four parameters are retrieved 
for each AIRS FOV: up to two layers each for both Tq and 
EOF. By retrieving each AIRS FOV individually, better fits 
between the simulated and observed cloud configurations are 
obtained. As a result, this adds an additional spatial resolu- 
tion of a factor of 3, and a sampling factor of 9, for Tq and 
leads to a higher degree of realism, especially for variability 
within complex cloud scenes and along cloud edges. If only 
one layer is retrieved, it is reported in the upper cloud layer, 
regardless of the altitude. 

The v6 retrieved Tq and ECF also benefit from other algo- 
rithm improvements in the cloud retrieval and elsewhere in 
the retrieval system. Since we are allowing for two horizon- 
tal layers of clouds as seen from above, there is the potential 
to create a nonexistent cloud layer that in practice only fits 
noise, compensates for other retrieval errors, or gives an un- 
physical best-fit solution of <0 or > 100% cloudiness. Our 
conjecture is that this “noise” in the ECF retrievals arises 
from the assumption of horizontally uniform surface temper- 
ature and emissivity fields, and the assumption of uniform 
temperature and water vapor profiles in the AIRS/ AMSU 
FOR. Any noncloud horizontal variability may be interpreted 
as additional structure in the ECF fields. The v6 algorithm 
better constrains the compensating cloud layer cases by im- 
proved logic for deciding when a one-cloud-layer solution 
adequately matches the radiances, and for converging to a so- 
lution away from a shallow minimum where the best-fit solu- 
tion is unphysical. Another improvement reduces the chance 
of placing the cloud above an inversion, which could prevent 
convergence to a solution. 

AIRS channels sensitive to low clouds are also sensitive 
to the surface. Thus, in overcast conditions, the surface tem- 
perature and emissivity are difficult to determine. In v5, we 
addressed this by using a microwave-only (AMSU) solution 
in about 17 % of cases to determine the clouds. In v6, the land 
surface emissivity is retrieved starting from a surface clima- 
tology derived from the monthly MODIS MYD11C3 emis- 
sivity product, and modified to fit the AIRS spectral channels 
using the MODIS baseline-fit emissivity approach (Seemann 
et al., 2008). This provides a more reasonable first guess 
and more stable solution when there is little information 
available. In v5, shortwave and longwave window channels 
were used simultaneously to retrieve the surface parameters 
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Table 1. Mean v5 ( v6 ) global (land and ocean only, and land+ocean) area- weighted cloud properties for January 2007. Cloud frequency is 
defined as an AIRS FOV that contains a value of ECF>0.01 (the sum of both cloud layers). 



Upper 

Lower 

Upper 

Lower 

Total 

Cloud 


Tq (K) 

Tq (K) 

ECF 

ECF 

ECF 

Frequency 

Global 

241.9 

271.5 

0.246 

0.205 

0.451 

80.5 % 

v5 ( v6) 

(243.4) 

(266.1) 

(0.282) 

(0.156) 

(0.438) 

(80.1%) 

Land 

235.9 

266.7 

0.210 

0.198 

0.408 

70.8% 

v5 (v6) 

(236.8) 

(265. 7) 

(0.240) 

(0.180) 

(0.420) 

(70.1%) 

Ocean 

244.9 

274.0 

0.264 

0.208 

0.472 

85.5% 

v5 (v6) 

(246.8) 

(266.3) 

(0.303) 

(0.145) 

(0.448) 

(85.4%) 


(7’ s and emissivity) often resulting in unstable solutions in 
the presence of clouds and highly reflective surface features 
(Hulley et al., 2009; Hulley and Hook, 2012). However, in v6 
only shortwave window channels are used to retrieve surface 
temperature (7’s), which results in more accurate determina- 
tion of spectral emissivity under more difficult cloud con- 
ditions (Susskind and Blaisdell, 2008). Furthermore, v6 uses 
an IR-microwave neural net solution (Blackwell et al., 2008) 
as the first guess for temperature and water vapor profiles 
and 7’s, which allows for reasonable solutions for many more 
cases than in v5. In the most overcast conditions over ocean 
(about 7 % of cases), v6 uses the neural net surface temper- 
ature directly when calculating clouds, resulting in a much 
better depiction of low clouds. 

Cloud retrievals for a single AIRS granule in the subtropi- 
cal western Pacific Ocean region on 6 September 2002 are 
shown in Fig. 1. This scene was selected as a representa- 
tive example because of the very large mix of cloud types 
and weather regimes found within it. The major weather fea- 
tures include tropical cyclone Sinlaku at the western edge, 
and a frontal system to Sinlaku’s north that separates a re- 
gion of broken and shallow cumulus, thin cirrus, and mul- 
tilayered cloud structures to the south, and more uniform 
stratocumulus-type layers in the north. The benefits of the 
single-FOV v6 cloud retrievals are apparent in the Tq struc- 
ture, especially along the edges of the frontal system and 
tropical cyclone, and within the cloud cover to the south of 
the frontal system. The frequency of FOVs containing two- 
layered clouds is significantly reduced in v6 compared to 
v5 (not shown), indicating an improvement in v6 by reduc- 
ing the frequency of compensating cloud layer cases. A few 
FOVs within the outer rain bands of Sinlaku contain unreal- 
istically warm Tq and low ECF, but the eye is much sharper 
than in v5 and has a more realistic diameter. 

The global mean properties of Tq, ECF, and cloud fre- 
quency are summarized in Table 1. At the AIRS FOV scale, 
about 80 % of the area of the globe at any given time contains 
cloud with ECF >0.01, and this value is about 70 % over land 
and 85 % over ocean. Most of the cloud amount is contained 
in the upper layer, and ocean FOVs contain more than twice 


the ECF in the upper layer (0.303) compared to the lower 
layer (0.145). The differences between v5 and v6 are rela- 
tively minor except for the increase (decrease) of ECF in the 
upper (lower) layer in v6, and a compensating change of a 
few K in Tq. The radiative consistency (Nasiri et al., 2011) 
of the cloud and surface products is nearly identical between 
v5 and v6 (not shown) and further implies the presence of 
compensating (and improved) changes in the ECF and Tq 
fields. 

Global distributions for both layers of Tq and ECF are 
shown for 2007 in Fig. 2. The tropical convective regions 
that contain cold cirrus are clearly depicted. The upper layer 
of Tq in the subtropical stratocumulus regions is significantly 
warmer than in v5 and indicates a substantial improvement in 
height assignment (to be quantified below). Furthermore, the 
lower layer of Tq is warmer in the low latitudes and sug- 
gests that improvements in v6 could lead to more realistic 
multilayered cloud configurations. This topic warrants fur- 
ther investigation. The majority of the cloud amount is con- 
tained within the upper layer as shown by the magnitudes of 
ECF. Interestingly, the lower layer of ECF is proportionally 
higher within tropical convection over South America, cen- 
tral Africa and the Intertropical Convergence Zone (ITCZ) 
compared to the Maritime Continent. Significant values for 
the lower layer ECF are also found in the storm tracks that 
are associated with nimbostratus clouds. This phenomenon 
also occurs in v5 and is a result of the tenuous nature of the 
upper portions of nimbostratus clouds and the tendency of 
the cloud retrieval algorithm to fit a second layer with a large 
value of ECF well within the cloud layer (Kahn et al., 2008a). 

2.2 AIRS, CloudSat, and CALIOP cloud top height 
histograms 

Figure 3 shows histograms of AIRS Zq over global 
oceans for v5 and v6 compared to cloud tops obtained 
from the CloudSat 94 GHz radar (Stephens et al., 2008) 
and the Cloud-Aerosol Lldar with Orthogonal Polarization 
(CALIOP; Winker et al., 2010) instruments for both sin- 
gle and two-layered clouds. The CloudSat cloud tops are 
taken from the release 4 (R04) 2B-GEOPROF product, 
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Fig. 1 . Granule maps from 6 September 2002. Upper row (left to right) is the 7b 1231 (K), the v6 upper-layer Tq, and the v5 upper-layer Tq. 
Middle row (left to right) is the v6 lower-layer Tq, the v6 total ECF, and the v6 cloud thermodynamic phase. Lower row (left to right) is the 
ice cloud optical thickness r, ice cloud effective diameter D e (pm), and ice cloud top temperature Tq ice (K). The granule number is 044 
with a start (end) time of 04:23:26 (04:29:26) UTC. 


and the CALIOP cloud top from the 5 km feature mask 
CAL_LID_L2_05kmCLay-Prov-V3-01 . The peaks in low- 
cloud top frequency observed by CloudSat and CALIOP 
are similar and at about 1.0-1. 5 km, with some slight di- 
urnal variability. In the AIRS single-layered retrievals (re- 
stricted to clouds with ECF >0.01), a broad peak in v5 
is located at 2-5 km, and in v6 lowers to 1.5-3. 5 km. In 
the AIRS v5 two-layered retrievals, the peak is located 
even higher than in the single-layered retrievals. However, 
in v6, the location of the peak is very similar between 
single- and two-layered retrievals. The diurnal differences 
in AIRS are minor for this subset of retrievals, although a 


slightly higher (lower) frequency of two-layered clouds are 
observed at night (day) in v6. There is a more prominent 
peak near the tropical tropopause in v6 if all values of ECF 
are included (not shown). Thus, by filtering out clouds with 
ECF <0.01, a much more realistic vertical distribution is ob- 
tained. CALIOP detects much more frequent tropical thin 
cirrus, although with a much wider peak than AIRS, and 
even more at night because of increased signal-to-noise ra- 
tio in CALIOP (Sassen et al., 2008; Winker et al., 2010). 
A small peak near 8 km in the single-layered v6 case and a 
broad peak from 6 to 10 km in the two-layered v6 clouds are 
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Fig. 2. The AIRS standard L2 cloud top temperature Tq (upper row) and effective cloud fraction ECF (lower row) for the upper layer (left 
column) and lower layer (right column) for 2007. 


Ocean Single-Layer V5 Ocean Single-Layer V6.0.2 




Land Single-Layer V5 Land Single-Layer V6.0.2 



Land Double-Layer V5 Land Double-Layer V6.0.2 



Counts Counts 


Fig. 3. Cloud top height Zq frequency for AIRS v5 (left column) 
and v6 (right column) over the ocean for single-layered cases (up- 
per row) and two-layered cases (lower row) for a one-week period 
in January 2007. The CloudSat and CALIOP collocated Zq is also 
shown on each panel. The CloudSat and CALIOP observations are 
very similar, but not exactly equal to each other, between the dif- 
ferent panels because of slight changes in the AIRS cloud detection 
between v5 and v6. 


seen in Fig. 3, and both are largely consistent with the active 
sensors. 

Figure 4 is similar to Fig. 3 except that the global land 
results are shown (only AIRS clouds with a two-layer 


Fig. 4. Same as Fig. 3, except for land. 

ECF >0.01 are included). Improvement for the low clouds 
over land is even greater than over the ocean. However, a 
sharp peak in high thin cloud is observed near 17 km unlike 
the ocean cases. An improvement in the realism of the broad 
peaks near 8-12 km is observed in v6 and is more consis- 
tent with the active sensors. In summary, the AIRS v6 cloud 
top height distributions are significantly improved over v5, 
although AIRS still struggles with locating the height of the 
thinnest cloud near the tropopause. 
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3 AIRS cloud thermodynamic phase and ice cloud 
property retrievals 

A new set of cloud products are added to the v6 support prod- 
uct files: these include cloud thermodynamic phase, ice cloud 
optical thickness (r), ice cloud effective diameter ( D e ), and 
ice cloud top temperature (Tcjce)- Jin (2012) and Jin and 
Nasiri (2013) describe the cloud thermodynamic phase algo- 
rithm, the theoretical basis for the channel selection, and the 
comparison of cloud thermodynamic phase between AIRS 
and CALIOP. The radiative transfer simulations of the ice 
cloud 7b spectra using the delta-four-stream (D4S) radiative 
transfer model (Yue and Liou, 2009) are described in Ou et 
al. (2013), while the optimal estimation retrieval approach is 
adapted from Bowman et al. (2006). 

3.1 Cloud thermodynamic phase 

The cloud thermodynamic phase is based on a series of T\, 
thresholds and A 7 j, tests that distinguish the spectral sig- 
natures of liquid and ice clouds (Nasiri and Kahn, 2008; 
Jin, 2012; Jin and Nasiri, 2013). These tests are applied 
to AIRS FOVs with a total two-layer ECF>0.01. A brief 
description of the key elements of the algorithm follows. 
There are four different ice-phase tests used: (1) 7j, at 
960 cm -1 (7^960) <235 K, which tests for the presence of 
optically thick and cold cirrus; (2) the 7b difference between 
1231 and 960 cm -1 (A7b,i23i-96o) >0.0K; (3) the 7b differ- 
ence between 1231 and 930cm -1 (A7b,i23i-96o) >1.75K; 
and (4) the 7b difference between 1227 and 960 cm -1 
( A7b, 1231— 96o) > — 0.5 K. Generally speaking, the colder and 
convective types of ice clouds pass more tests, and thin cirrus 
passes fewer tests (discussed later). Although detection of ice 
within the AIRS FOV is the main purpose of the algorithm, 
there are two different liquid-phase tests: (1) the 7b dif- 
ference between 1231 and 960cm -1 (A7b.i23i-96o) < — 1.0 
K, and (2) the 7b difference between 1231 and 930 cm -1 
(A7bq23i-96o) < — 0.6 K. All of the individual test results 
are reported in the field cloud _phase_bits. Similar to the 
ice cloud cases, the uniform and homogeneous liquid clouds 
more often pass both tests than do broken clouds. In the case 
of “unknown” cloud, Jin and Nasiri (2013) show that over 
99% of these cases pass no ice or liquid tests (less than 1 % 
have at least one liquid and ice test each simultaneously pass- 
ing). After summing the results of all tests, “ice” is obtained 
if the value is positive (between +1 and +4) and "liquid” if 
negative (either —1 or —2). These values are reported in the 
field cloud _phase_3x 3. 

The granule map of cloud thermodynamic phase in Fig. 1 
shows numerous ice clouds with phase values from +1 to 
+4. The frontal band to the north contains values from +1 
to +3, while a value of +4 shows up in small areas in the 
rain bands around tropical cyclone Sinlaku. A significant 
amount of unknown phase (0) is found to the east of Sinlaku 
and south of the frontal band. These unknown-phase clouds 


are low in altitude (high Tq), are most likely broken given 
the low values of ECF retrieved, and most closely resemble 
shallow trade cumulus: marine boundary layer clouds of this 
type most often populate the unknown-phase category. Liq- 
uid clouds with values of — 1 and —2 show up north of Sin- 
laku and the frontal band. A cloud-phase value of —2 is asso- 
ciated with higher values of ECF, and shows that a stronger 
spectral signature is obtained from more homogeneous and 
optically thicker clouds (e.g., Nasiri and Kahn, 2008; Kahn 
et al., 2011). 

Figure 5 shows global patterns of cloud thermodynamic 
phase for 2007. Ice cloud frequencies approach or exceed 
80-90% over the Maritime Continent and adjacent regions, 
the ITCZ, central Africa, and the tropical portions of South 
America. These climatological patterns are most similar to 
those previously obtained from HIRS (Wylie et al., 2005; 
Stubenrauch et al., 2006), AIRS (Stubenrauch et al., 2010), 
and other passive sounders (Liou, 1986), with higher magni- 
tudes compared to CALIOP (Hu et al., 2009), MOD1S (Yang 
et al., 2007; Meyer et al., 2007), and 1SCCP (Rossow and 
Schiffer, 1999). These differences in frequency are mostly 
explained by the data-filtering approach (colder than —40 °C 
in Sassen et al., 2008), the pixel size (larger pixels are more 
likely to contain an ice signature), and instrument or algo- 
rithm sensitivities (e.g., Ackerman et al., 2008). 

Despite the hemispheric differences in ice and liquid cloud 
frequencies, the overall cloud frequency (sum of liquid, ice, 
and unknown) is very similar between the two oceanic storm 
tracks (Fig. 5). These values are much larger than the ECF 
values shown in Fig. 2 because the phase algorithm only tests 
for the presence of cloud within the AIRS FOV rather than 
its fractional area or opacity. As expected, there is a high 
frequency (50-60%) of liquid clouds in the subtropical stra- 
tocumulus regimes near the continents. However, in compar- 
ison to ISCCP, MODIS, and other data sets, the magnitudes 
of liquid frequency are lower and, correspondingly, the mag- 
nitudes of unknown frequency are higher. This is not surpris- 
ing given that the phase retrieval depends only on an infrared 
spectral signature. Future algorithm improvements may re- 
duce the fraction of unknown cloud categorization. 

Stratocumulus clouds cover approximately 20 % of Earth’s 
surface, but coverage can be as high as 25-40 % over mid- 
latitude oceans (Wood, 2012). AIRS’ limitations in low- 
cloud categorization are very similar to those from HIRS 
(e.g., Wylie et al., 2005). Approximately 60% of all liquid 
clouds are identified by AIRS as unknown when compared 
to CALIOP cloud phase (Jin, 2012; Jin and Nasiri, 2013). 
Visual inspection shows this value is higher for trade cumu- 
lus and lower for stratocumulus. In general, the AIRS phase 
product is very conservative in detecting liquid clouds; rarely 
is an ice cloud (according to CALIOP) identified as liquid by 
AIRS (Jin, 2012). The AIRS IR phase tests are much more 
sensitive to ice compared to liquid, partly because of the ther- 
mal surface contrast for high cold clouds (e.g., Kahn et al., 
2011). A high frequency (50-70 %) of liquid water clouds is 


www.atmos-chem-phys.net/14/399/2014/ 


Atmos. Chem. Phys., 14, 399-426, 2014 


406 


B. H. Kahn et al.: The Atmospheric Infrared Sounder version 6 cloud products 




Unknown Froquoncy 



Liquid Froquoncy 



Fig. 5. Cloud thermodynamic phase for ice (upper left), liquid (lower left), unknown phase (upper right), and the sum of the three phases 
(lower right) for 2007. The relative global frequencies of occurrence (% of total AIRS FOVs) are 26.5 % for ice, 16.2 % for liquid, and 34.2% 
for unknown. 


observed in the austral summertime between 50° S and the 
Antarctic coast (Wood, 2012). These high frequencies are 
also consistent with supercooled liquid cloud frequency ob- 
served by the Lidar In-space Technology Experiment (LITE; 
Hogan et al., 2004) and CALIOP (Hu et al., 2010; Choi et al., 
2010). Similar patterns around Antarctica are also observed 
in other austral summertime months and years, and liquid 
water frequencies are drastically reduced during austral win- 
ter (not shown). 

The spatial distributions of the frequency of passed liquid 
and ice tests are shown in Fig. 6. Values of —2 (both liq- 
uid cloud tests passed) show up most frequently along the 
Antarctic coast, within the confines of the coastal subtropi- 
cal stratocumulus regimes, and in a small portion of south- 
east Asia (e.g., Wood, 2012). Values of — 1 (one liquid cloud 
test passed) are more prominent in the subtropics between 
the two storm tracks, and also westward of the coastal sub- 
tropical stratocumulus regimes. A weaker spectral signature 
implies a lower likelihood that both liquid tests pass, and val- 
ues of — 1 correspond closely to the presence of broken liquid 
clouds. Positive values indicate passed ice cloud tests. Values 
of +1 are most prominent over land areas and the tropical 
Pacific and Indian oceans. Values of +2 are found in simi- 
lar regions with some additional emphasis on the midlatitude 
storm tracks. Values of +3 are found in a more confined area 
of the tropics and smaller areas of the midlatitudes. Values of 
+4 are the least common and show up within deep convec- 
tion in the tropical western Pacific, tropical South America, 
and central Africa, with negligible frequencies elsewhere. 


3.2 Ice cloud property retrievals 

The new ice cloud property retrievals r, D e , and 7c, ice are 
shown in Fig. 1. The 7c , ice and upper-layer Tq correspond 
well with each other, although warmer values are observed 
in places where the AIRS retrieval obtains two layers. This 
suggests that ice cloud retrievals may be improved in the fu- 
ture by including an additional layer(s). The optically thick- 
est cirrus is located near the center of Sinlaku and in the core 
of the frontal band. The D e pattern is more variable (imply- 
ing higher noise) than r and 7c, ice and a higher frequency of 
data dropouts exists because of the more stringent QC com- 
pared to r and 7c, ice- However, many distinct granule-scale 
patterns emerge. An abundance of cirrus to the south of Sin- 
laku is observed to have much lower r compared to the ice 
clouds surrounding the eye. 

3.2.1 Ice cloud property retrievals 

The three primary ice cloud products are retrieved simul- 
taneously and are found in the AIRS L2 support prod- 
uct files: r ( ice_cld_opt_dpth ), D e ( ice_cld_eff_diam ), and 
7c, ice {ice _cld temp _eff). 7c, ice is considered a less crucial 
parameter than r and D e since it is already retrieved as a part 
of the AIRS L2 standard product suite, but much improved 
X 2 fits and more frequent convergence were obtained when 
retrieving 7c ice as a third parameter. The first two parame- 
ters (r and D e ) are retrieved in log space to prevent negative 
values. The three parameters are simultaneously retrieved 
on individual AIRS FOVs that contain ice (anywhere from 
+ 1 to +4 tests passed) using an optimal estimation retrieval 
as a post-processor after the AIRS standard L2 retrieval is 
completed. The optimal estimation algorithm is derived from 
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Fig. 6. The various phase tests for liquid (upper row) and ice (middle and lower rows) phases for 2007. In the case of liquid, —2 indicates 
that the two liquid-phase tests were passed, and —1 indicates only one of the two tests passed. In the case of ice phase, maps are shown for 
the number of ice-phase tests that passed (1-4). The relative frequencies of occurrence (% of total AIRS FOVs) for tests —2, —1, 1, 2, 3, and 
4 are 7.9 %, 8.3 %, 6.7 %, 10.7 %, 8.2 %, and 0.8 %, respectively. Please refer to Sect. 3. 1 for an explanation of each phase test. 


the Tropospheric Emission Spectrometer (TES) retrieval, de- 
scribed in Bowman et al. (2006). The algorithm minimizes 
the cost function 

C=\\y-F(x,b)\\ 2 ,+Wx-XaW 2 ,, (1) 

where y is the vector of measured radiances, Fix) is the 
forward-modeled radiance vector, jr is the state vector of the 
parameters retrieved (r, D e , and 7’c , [(_. K ) , x a is the a pri- 
ori state vector, b is a vector containing fixed atmospheric 
state variables (temperature profile, water vapor profile, etc.) 
and observational metadata necessary for calculating the ra- 
diances, S” 1 is the inverse noise covariance (diagonal for 
AIRS-footprint retrievals), and S“ 1 is the inverse covariance 
of the a priori. 

For this initial effort, construction of the x a and S“* ma- 
trices is mostly ad hoc. The inputs to our a priori state vec- 
tor are an assumed r = 3.0 and D e — 30 pm, while 7c, ice is 
initially from the AIRS L2 standard product upper-level 7/ . 
The fixed first-guess values for r and D e were settled on af- 
ter investigating more dynamic initial guesses that depended 
on 7c and other parameters. These other approaches led to 
poorer spectral radiance fits and reduced occurrences of re- 
trieval convergence. 


S“* contains the a priori variances along the diag- 
onal, with zero off-diagonals, and are also assumed to 
be constants ( log_ice _cld_opt_dpth _prior_var = 0.111, 
log ice cld eff diam _prior_yar = 0.16, and ice_cld_ 
temp _ejf _prior _var — 225). These a priori values are consis- 
tent with reported histograms of these cloud properties from 
remote sensing observations. The variances are in practice 
dependent on many physical variables including cloud and 
scene type, season, latitude, altitude, and possibly many other 
factors. Whether independent satellite, in situ, or model- 
derived data sets are the appropriate proxy for prior first 
guesses and variances is uncertain and warrants further re- 
search. Thus, with this use of assumed constants in the a pri- 
ori, we caution against quantitatively using the reported error 
or averaging kernels (described below) until further research 
is performed. 

The solution of the above cost function follows the Gauss- 
Newton and Levenberg-Marquardt methods described in 
Section IV of Bowman et al. (2006). The line-by-line ra- 
diative transfer model (LBLRTM) forward model was re- 
placed with the stand-alone AIRS radiative transfer algo- 
rithm (SARTA; Strow et al., 2006), and the cloudy sky spec- 
tra are calculated from a version of SARTA that is coupled to 
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a delta-four-stream (D4S) cloudy radiative transfer program 
to account for ice cloud scattering (Ou et al., 2013). 

Ignoring errors from fixed atmospheric state variables, and 
assuming that the retrieval is somewhat linear in the neigh- 
borhood of the solution (see Sect. 5.4 in Rodgers, 2000), we 
calculate the solution error covariance as 

S= (k^-'K+S- 1 ) , (2) 

where K is the Jacobian (dy/djr), calculated by finite dif- 
ferences. We assume the errors are uncorrelated, and the re- 
ported errors are the square roots of the diagonals of S. Note 
that the retrievals for r and D e are perfonned in natural log 
space but are reported in the output files ( ice_cld_opt_dpth 
and ice_cld_eff_diam) in linear space. However, the reported 
errors ice_cld_opt_dpth_err and ice_cld_eff_diam_err have 
been left in natural log space; that is, 

ice_cld_opt_dpth_err = s (log e [r]) (3) 

and 


ice_cld_eff_diam_err — e (log e [7) e ]) ■ (4) 

The lower (upper) boundaries for the retrieved r and D e 
can be determined by dividing (multiplying) by the exponen- 
tial of the reported error fields. For example, the range [lower 
error, higher error] of the retrieved r can be calculated by 

{e -ice_cld_opt_dpth_err x ice _ c ld_ opt _d P th, (5) 

e ice_cld_opt_dp,h_err x ice _ c ld_op t _dpth], 


and the same approach is used for D e . The reported error for 
7c , ice ( ice cld temp eff err ) is in linear space. Again, for 
reasons given, we caution against using these errors quanti- 
tatively. Further discussion regarding error characterization 
is addressed in Sect. 3.2.4. 

The averaging kernel matrix, A, is a measure of the sensi- 
tivity of the retrieval to a change in the true state: 

A = ^ retrieved = (Vs-'K + S^rVs-'K. (6) 

The reported scalar averaging kernels (AK) 
ice_cld_opt_dpth_ave_kern , ice_cld_eff_diam_ave_kern, 
and ice_cld_temp_eff_ave_kem are from the diagonal of 
A, and independently range from 0.0 to 1.0. The / 2 fitting 
parameter is calculated as 


X 


2 


J_y / yj-[F(x)]j 
N 2—ii=l ^ g. 


(7) 


where e,- is the radiance error in channel i and N is the num- 
ber of channels. 


3.2.2 QC and channel selection 

The three most important ancillary parameters are the qual- 
ity control (QC) indicators for r ( ice_cld_opt_dpth_QC ), D e 
{ice_cld_eff_diam_QC), and 7c, ice {ice _cld temp _eff_QC). 
For ice_cld_opt_dpth_QC and ice_cld_temp_eff_QC ? the 
range is from 0 to 2, where 0 = “best”, 1 = “good”, and 
2 = “do not use”. We strongly discourage the use of scenes 
with QC = 2 unless users carefully validate the retrieval re- 
sults or consult with members of the AIRS Science Team. For 
ice_cld_eff_diam_QC, only values of 1 and 2 are reported. 
Since D e is a very challenging parameter to retrieve and in- 
terpret, a conservative QC approach was decided to be most 
appropriate. The r, D e , and 7c, ice parameters must be used 
in conjunction with the QC at the FOV scale, since values 
with QC = 2 are also reported in the AIRS L2 support prod- 
uct. These products are also available as gridded level 3 (L3) 
files with appropriate QC applied. 

The QC indicators are derived from the 
ice_cld _fit_reduced_chisq (x 2 ) between the simulated and 
observed 7b’s and scalar AKs ice_cld_opt_dpth_ave_kern, 
ice_cld_eff_diam_ave_kern, and ice_cld_temp_eff_ave_kern 
that independently range from 0.0 to 1.0. The QC indicators 
derived from combinations of x 2 mid AKs for the three 
retrieval parameters are described in Table 2, and additional 
retrieval parameters are listed in Table 3. All values are 
reported in the AIRS L2 support product files. These QC 
indicators are neither absolute nor quantitative, but should be 
interpreted as an approximate indicator for robust retrieval 
values that are obtained from good spectral 7b fits between 
simulations and observations and higher values of informa- 
tion content. As the QC indicators for these fields are not 
identical, it is possible that one parameter may have QC = 
0, while the other parameters may have QC = 1 or 2. The 
percentage occurrence of each QC value for each retrieval 
parameter during 2007, where 27.0 % of all AIRS FOVs are 
identified as containing ice, is shown in Table 4. 

An illustrative retrieval for a thin cirrus cloud (r = 0.46, 
77 e = 41(im, and 7c ice = 213 K) located at 15.6° N and 
132.6° E is shown in Fig. 7. There are 59 AIRS channels in 
the 8-15 pm window region that are used for the retrieval 
(see Appendix C), and are different from those used in Kahn 
et al. (2008b). The channels were not optimized for maximiz- 
ing the information content of the retrieved cloud parameters 
nor were the IR spectra compressed to reduce noise (Huang 
and Antonelli, 2001; Antonelli et al., 2004) or to increase 
computational efficiency of the radiative transfer (Liu et al., 
2005). Strong water vapor absorption lines were avoided, and 
channels with relatively low values of NEdT were retained. 
A smaller set of channels was tried and resulted in a reduced 
frequency of retrieval convergence. Likewise, for a larger set 
of channels, the computational expense increased beyond ac- 
ceptable values. The sensitivity to adjustments in r, D e , and 
7c, ICE is also shown in Fig. 7, and as previous investigations 
have revealed, the greatest effect on the 7], spectrum is due 
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Table 2. Listed is the quality control (QC) determination for the three ice cloud retrieval parameters. The scalar averaging kernels (AK) are 
for each parameter and the value of x 2 is the measure of fit between the observed and simulated radiances. The QC indicators are neither 
absolute nor quantitative. The variables ice_cld_ftt_reduced_chisq (x 2 ), icecldoptdpthavekern (AK for t), ice_cld_eff_diam_ave_kern 
(AK for D e ), and ice_cld_temp_eff_ave_kern (AK for 7c,ICe) are reported in the AIRS L2 support product files. 



Best 

(QC = 0) 

Good 
(QC = 1) 

Do Not Use 
(QC = 2) 

Ice Cloud Optical 
Thickness (r) 

Both AK> 0.8 
and x 2 < 10.0 

Either AK> 0.8 
or x 2 < 10.0 

Both AK <0.8 
and x 2> 10.0 

Ice Cloud Effective 
Diameter ( D e ) 


Both AK >0.8 
and x 2< 10.0 

Either AK >0.8 
or x 2 < 10.0 

Ice Cloud Top 
Temperature (7c,ICe) 

Both AK >0.8 
and x 2 < 10.0 

Either AK>0.8 
or x 2 < 10.0 

Both AK <0.8 
and x 2> 10.0 


to changes in r. For a fixed value of r, the Icjce is not very 
sensitive to height changes, but the highest sensitivity is seen 
in the C02-slicing channels. With x and 7 cjce fixed, D e 
shows the least sensitivity of the three parameters. The well- 
known sensitivity in slope across the atmospheric windows is 
best observed in this example when comparing D e = 20 pm 
with the other values of D e . This subtle sensitivity shows the 
challenge of retrieving D e , especially for clouds with low (or 
high) values of r. This demonstrates that a robust retrieval 
methodology like that presented in Bowman et al. (2006) is 
preferable to an ad hoc look-up table approach (Kahn et al., 
2008b). 

Some evidence of the ice_cld_opt_dpth _Jirst_guess re- 
mains in the results to follow, although only a negligible per- 
centage of retrievals stick to a value near the first guess. For 
the ice_cld_eff_diam _first_guess, there is no evidence that 
the retrieval sticks to near the first guess value. However, 
this parameter is sensitive to the width of the finite differ- 
ence that is fixed at 10 pm, which is the D e interval of the 
single-scattering properties (Baum et al., 2007). When much 
finer binning is performed in the plotting of D e , for instance 
at 1 pm intervals (not shown), a much higher frequency of oc- 
currence near the 10 pm increments is found than in between 
the increments. 

3.2.3 Global distributions 



Wave Number (cm ) 


Fig. 7. Shown is an example AIRS spectrum on 6 September 2002 
located in granule #44 (same as Fig. 1). The retrieval is for a rel- 
atively thin cirrus cloud and is located at 15.6° N and 132.6° W. 
The upper panel shows the best fit for r = 0.46, and adjustments 
to 1.0, 2.0, 3.0, and 4.0. The middle panel shows the best fit for 
D e = 41 pm, and adjustments to 20.0, 60.0, 100.0, and 140.0 pm. 
The lower panel shows the best fit for Tq jqj = 213 K, and adjust- 
ments to 200, 220, and 240 K. 


The global ice cloud properties for 2007 are shown in Fig. 8. 
The r distributions in the tropical western Pacific Ocean con- 
tain a narrow band of high values compared to ice cloud 
frequency. The highest values of r are associated with the 
convective band closest to the 1TCZ, while this region is sur- 
rounded by a wide latitudinal extent of cirrus with lower val- 
ues of r. Furthermore, r is higher in the winter oceanic storm 
tracks (Sect. 3.2.6), and is greatly reduced in the Arctic re- 
gion (Curry and Ebert, 1992) and also over the cold win- 
tertime east Asian and North American continental regions. 
Very low values of r dominate the subtropical subsidence re- 
gions. 


The D e distributions are complex and appear to be non- 
intuitive. They have a broad minimum in the high latitudes, 
corresponding closely to r. D e also has a minimum over the 
tropical western and central Pacific Ocean, especially on ei- 
ther side of the 1TCZ, where thin cirrus is most common. 
A small maximum in D e along the 1TCZ is consistent with 
larger MODlS-derived r c obtained within deep convective 
events in the tropics (Yuan and Li, 2010; cf. Fig. 12). Re- 
trievals of r e from surface-based radar and lidar observations 
also show somewhat larger values in profiles classified as 
more convectively active (Protat et al., 2011). These results 


www.atmos-chem-phys.net/14/399/2014/ 


Atmos. Chem. Phys., 14, 399-426, 2014 


410 


B. H. Kahn et al.: The Atmospheric Infrared Sounder version 6 cloud products 


Table 3. Names and descriptions of all ice cloud retrieval parameters found in the AIRS L2 support product files. 


AIRS L2 support product variable 
name 

Description 


ice cld opt dpth 

Ice cloud optical thickness (r) in the 
retrieval state vector x 


ice cld eff diam 

Ice cloud effective diameter (D e ) in the 
retrieval state vector x 


ice _cld temp eff 

Ice cloud top temperature (Tqice) in the 
retrieval state vector x 


ice cld opt dpth QC 

Quality control (QC) flag for r (see Table 2) 


ice cld eff diam QC 

QC flag for D e (see Table 2) 


ice cld temp eff QC 

QC flag for 7c, ICE ( see Table 2) 


log ice cld opt dpth _prior_var 

A priori variance for r from diagonal of ST 1 . 
fixed to constant value of 0. 1 1 1 


log ice cld eff diam _prior_yar 

A priori variance for D e from diagonal of 
ST 1 , fixed to constant value of 0. 1 6 


ice_cld_temp_eff_prior_var 

A priori variance for Tq ice from diagonal of 
ST 1 , fixed to constant value of 225 



icecldoptdpthavekern Scalar averaging kernel (AK) for r from 

diagonal of A, ranges from 0.0 to 1.0 


ice cld eff diam ave kern 

Scalar averaging kernel (AK) for D e from 
diagonal of A, ranges from 0.0 to 1.0 

ice _cld temp eff _ave kern 

Scalar averaging kernel (AK) for Tcjce from 
diagonal of A, ranges from 0.0 to 1.0 

ice cld opt dpth Jirst_guess 

First guess of r, fixed to constant value of 3.0 in the 
a priori state vector x a 

ice_cldjemp_effjirst_guess 

First guess of 7c, ICE' varies and is set to the upper- 
level AIRS L2 Tq in the a priori state vector x a 

ice cld eff diam Jirst_guess 

First guess of r, fixed to constant value of 30 pm in 
the a priori state vector x a 

ice cld opt dpth err 

Error estimate for r from diagonal of S 

ice cld eff diam err 

Error estimate for D e from diagonal of S 

ice cld temp eff err 

Error estimate for Tc.ICE from diagonal of S 

ice cld Jit reduced chisq (x 2 ) 

Chi-squared fit between observed and simulated 
AIRS radiances (see Eq. 4) 


Table 4. Percent occurrence of the QC values (0, 1 or 2) for r, D e , 
and 7c jce f° r 2007. These percentages are from the subset of AIRS 
FOVs that contain ice (26.5 % globally during this time period). 



Best 
(QC = 0) 

Good 
(QC = 1) 

Do Not Use 
(QC = 2) 

r 

64.1% 

29.0% 

6.9% 

»e 


68.1% 

31.9% 

7c,ICE 

75.5% 

20.0% 

4.5 % 


appear to contradict retrievals from HIRS (Stubenrauch et 
al., 2004; cf. Fig. 7). However, the discrepancies with HIRS 
are probably a result of sampling, binning, and other differ- 
ences between this study and Stubenrauch et al. (2004). FOV- 
scale comparisons of retrievals between the NOAA HIRS 
and Aqua AIRS instruments are likely to reveal more mean- 
ingful information. Maximum values of D e are obtained on 
the subtropical side of the midlatitude jets in both hemi- 
spheres over the oceans. There is a decrease in D e poleward 
across both storm tracks. There is also some structure in ice 
cloud frequency and D e associated with orography in coastal 
Antarctica, the Andes, and other regions of elevated topog- 
raphy. A much larger data set spanning many years and all 


Atmos. Chem. Phys., 14, 399-426, 2014 


www.atmos-chem-phys.net/14/399/2014/ 


B. H. Kahn et al.: The Atmospheric Infrared Sounder version 6 cloud products 


411 


lc* Old Fr*q lc« Cld r 



Fig. 8. Ice-phase frequency (upper left, repeat from Fig. 5), r, D e (pm), and 7c, ICE (K) for 2007. 





Fig. 9. Histograms of r (left), 7c , ice (middle), and D e (right) for 2007. The FOV-scale retrievals that passed 1—4 ice tests are shown 
separately, as well as histograms for all retrievals (sum of tests 1-4). 


seasons is necessary to quantify the behavior of orographic 
clouds and relate their properties to static stability, vertical 
velocity, and flow regime (Joos et ah, 2010). 

The 7c, ice pattern is similar to the upper-layer 7) in 
Fig. 2, but is on average colder. This is expected because 
7c, ice is a subset of only ice clouds from the ensemble of 
all cloud types reported in 7c. 

3.2.4 Histograms sorted by cloud thermodynamic phase 
flag 

Histograms of the AIRS ice cloud parameters are shown in 
Fig. 9. The r histogram for all ice tests shows a broad dis- 
tribution with a peak occurrence near 0.5-0. 6. The a priori 
value of 3.0 shows up as a small bump. From this, we con- 
clude that the a priori guess for r has little impact on the 
broader retrieval results. For an ice-phase value of +1, the 
broad peak in r is located closer to 0.2-0. 3, and a reduction 
in occurrence frequency is seen for r > 1.0. The histogram for 
+2 has a peak closer to r = 0.4-0. 5 and a larger occurrence 
frequency is observed for r > 1.0 compared to +1. In the case 
of +3, there is a broad peak for r = 0. 7-1.0, and a smaller 
secondary peak for r = 6-7. For +4, the primary peak is 


near r = 5-6 and fewer values are observed for r < 1.0. In 
summary, at higher values of r, the spectral signature for ice 
is stronger, and more ice-phase tests are passed. Conversely, 
more tenuous ice clouds pass fewer tests, and are thus de- 
tected less robustly. 

A similar pattern is observed for 7c, ice- The ice clouds 
are colder on average as more ice tests are passed. There 
is a small bump in the occurrence frequency around 290 K 
for + 1 and +3 that suggests biases from multilayered cloud, 
compensating errors in the retrievals, and/or multiple solu- 
tions for very low r (e.g., Kalin et al., 2008b). The warn 
values of 7c, ice are most prevalent in the subtropics (Fig. 8). 
For ice tests +1 to +3, the peak is near 230 K, with a small 
peak near 2 1 5 K for + 1 . An additional subtle change in slope 
around 195 K for +1 and +2 indicates a signature of high al- 
titude and optically thin cirrus near the tropical tropopause 
layer (TTL). However, these cases do not dominate the over- 
all occurrence frequency. Given the inherent limitations in 
retrieving thin cirrus from thermal contrast observations, the 
single layer assumption, and uncertainties in the a priori at- 
mospheric state (e.g., Posselt et al., 2008), the AIRS retrieval 
is severely underestimating the occurrence of very thin TTL 
cirrus (Kalin et al., 2008b). CALIOP is better designed for 
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Fig. 10. Histograms of r (top row), 7c, ICE (middle row), and D e (bottom row) for QC = 0 (“best”, left column), QC = 1 (“good”, middle 
column) and QC = 2 (“do not use”, right column) retrievals for 2007. The histograms are organized by latitude band: tropics (30° S-30° N), 
SH (60-30° S) and NH (30-60° N) midlatitudes, and SH (90-60° S) and NH (60-90° N) polar regions. 


investigations involving TTL ice clouds (e.g., Sassen et al., 
2008, 2009). 

The histograms in Fig. 9 behave differently for D e com- 
pared to r and Tcjce in that their shapes are very similar 
among all ice test combinations. The peak frequency of D e 
occurrence varies from 40 to 50 pm and drops off substan- 
tially at smaller and larger diameters, and a gap in the re- 
trievals around 1 60 pm is apparent. Values in this size range 
are retrieved but are flagged as bad QC (see Fig. 10). The 
cause of this is uncertain and warrants further investigation. 

3.2.5 Histograms sorted by QC, latitude band, and 
error estimates 

The histograms of r, D e , and 7c, ice are sorted by QC and lat- 
itude band in Fig. 10. For QC = 0, the tropics are dominated 
by thin cirrus with a peak occurrence frequency of r = 0.2- 
0.3, and a much smaller but notable peak near r = 6.0, con- 
sistent with large particles lifted by convection. The polar 
areas lack high values of r as predicted by Curry and Ebert 
(1992). The midlatitudes are similar to the polar areas for QC 
= 0 but with slightly higher occurrences of larger r. For QC 
= 1 , there is a shift in the maximum occurrence frequency to 
lower values of r . However, the relative ordering of different 
latitude bands are similar for QC = 0 and QC = 1 . The sig- 
nature of the a priori (r = 3.0) is much more variable for QC 


= 1 (maximum in Antarctica, absent in tropics; we speculate 
this difference may be related to reduced information content 
in scenes with weaker thermal contrast), while the peak is 
small but consistent across all latitude bands for QC = 0. For 
QC = 2, the peak occurrence frequencies are at low values 
of r, which is sensible because these values are more suscep- 
tible to bad fits in scenes with multilayered clouds and low 
information content. For QC = 0 values of 7c, ice. the cold- 
est clouds are found in the tropics, and very warm clouds are 
found along the fringes of the tropics (see Fig. 6). The extra- 
tropics have a more confined 7c,ice distribution, consistent 
with a warmer tropopause and colder surface than the trop- 
ics. For QC = 1 and QC = 2 retrievals of Tcjce, the peak 
occurrence frequency is located at somewhat colder values 
in the extratropics, but is flatter and warmer in the tropics. 
The QC = 2 retrievals of Tcjce show several peaks and high 
counts of warm Tcjce, indicative of poor fitting and low in- 
formation content. QC = 1 and QC = 2 retrievals of D e peak 
between 40 and 50 pm for all regions. 

The fixed values for initial guesses, prior constraints and 
variances, and the absence of off-diagonal terms may ad- 
versely impact the magnitudes and dynamic ranges of the er- 
ror estimates calculated by Eq. (2). Despite these shortcom- 
ings in the v6 cloud retrieval algorithm, a qualitatively rea- 
sonable set of error estimates is obtained. These are presented 
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Fig. 11. Shown is the relative error (%) of r (left) and D e (right) as 
a function of r. Only retrievals with QC = 0 (r) and QC = 1 (r and 
D e ) are included. The color scale indicates the relative density of 
occurrences. 


as relative errors for r and D e and are shown in Fig. 11. The 
relative error for r decreases from 10 % to 2 % as the magni- 
tude of r increases from 0.1 to 1.0 and is somewhat constant 
for values of r > 1.0. This is consistent with the TES opti- 
cal depth error estimates in Eldering et al. (2008). There are 
small populations of retrievals with relative errors between 
20 and 50 % near r = 1 .0, and errors between 1 and 5 % 
for very thin cirrus with r <0.1. With regard to D e , most 
cloud retrievals have relative errors of approximately 10% 
for r = 0.1, reaching a minimum of 1-3 % near r = 1.0. As 
r increases further, the relative error increases to 5-10% for 
most values. Throughout the range of r, a small number of 
D e retrievals have relative errors greater than 10 %. Previ- 
ous sensitivity studies (e.g., Cooper et al., 2007; Kahn et al., 
2008b; Posselt et al., 2008) suggest much larger uncertain- 
ties on the order of 30-50% are expected for ice cloud D e . 
These previous studies included the impacts of uncertainties 
in atmospheric and surface state, ice crystal size and habit 
distribution, and the vertical geometry of clouds, among oth- 
ers, in retrieval uncertainty estimates of D e . None of these 
aforementioned factors are included in the present error es- 
timates. Therefore, the fact that the estimates presented in 
Fig. 11 are low is no surprise. 

3.2.6 Zonal averages 

Zonally averaged histograms of r, D e , and 7c, ice during 
2007 are shown in Fig. 12 and the seasonal variations are 
detailed in Fig. 13. All retrievals are sorted by land, ocean, 
day, and night. The oceanic r is highest in boreal winter near 
40° N and is greatly reduced at poleward latitudes. Another 
broad peak of r is found in the austral summer with a re- 
duction near Antarctica. However, the latitudinal r gradient 
is smaller in the summer SH; this pattern is also observed 
in other years (not shown). In the boreal wintertime midlat- 
itudes, r is slightly higher during day over ocean compared 
to land. The reverse is true in the austral summertime sub- 
tropics. There is a minimum of r near 10° N that is also ob- 
served in the MODIS Collection 5 (C5) December-January- 


February (DJF) time frame over both land and ocean (Hong 
et al., 2007). For AIRS, r over land is higher in the austral 
summer and may indicate an increased rate and/or vigor of 
convection, also in agreement with Hong et al. (2007). In 
Fig. 12, there is a pronounced minimum of r near 10-20° N 
over land. Day and night r values differ by as much as a fac- 
tor of 1.5, but the diurnal signal is smaller over ocean than 
over land. In the next section, we will show that the diurnal 
cycle needs additional regional context to fully describe the 
complexity and amplitude within the variety of cloud param- 
eters. At smaller scales in the presence of complex topogra- 
phy, the diurnal variations are much larger in magnitude than 
found in the zonal means. 

7c, ice has a very strong diurnal signal over land and less 
so over ocean, and reaches a maximum value in the subtrop- 
ics. This is also true in other years (not shown), although the 
magnitude varies from year to year. 

D e has a minimum in the tropics and polar winters, and 
a smaller minimum in Antarctica during summer. In the 
MODIS C5 data set, there is a maximum in r e during DJF 
over the deep tropical ocean and a weak minimum over land 
(Hong et al., 2007). In AIRS, there is a pronounced max- 
imum in the midlatitudes near 30^10° N/S on the equato- 
rial side of the jet stream, with a drop-off poleward of these 
latitudes. These results are most similar to the CAM5 con- 
trol run in Gettelman et al. (2010), while other CAM5 ex- 
periments tend to exaggerate the higher magnitude of r e in 
the midlatitudes compared to AIRS. Although Gettelman et 
al. (2010) report a cloud top value of r e , a more detailed com- 
parison taking into account the weighting functions of AIRS 
is warranted. The tendency for maxima in the midlatitudes 
and minima in the tropics and high latitudes are also seen 
in other climate model experiments. This includes simula- 
tions of transparent ice cloud in the 5th generation of the 
Max Planck Institute for Meteorology atmospheric general 
circulation model (ECHAM5) (Joos et al., 2010; cf. Fig. 5), 
and in modified versions of the CAM3 with various ice nu- 
cleation physics parameterizations (Wang and Penner, 2010; 
Liu et al., 2012). 


4 Additional tests 

Following are two process-based “stress tests” that are de- 
signed to gain further insight into the initial performance of 
the AIRS cloud products. The first test quantifies the varia- 
tions of key cloud parameters at two local times (01:30 and 
13:30 LT) over the Maritime Continent, where there is a very 
pronounced diurnal cycle in convective ice cloud and rainfall 
(Neale and Slingo, 2003; Nesbitt and Zipser, 2003; Tian et 
al., 2006; Qian, 2008). The second test composites midlati- 
tude cyclones in the two hemispheres using a previously pub- 
lished methodology by Naud et al. (2006, 2010). Substan- 
tial differences in cloud structure are found between the NH 
winter and SH summer storm tracks, which are exaggerated 
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Fig. 12. Zonal averages of r (left), 2c, ICE (center), and D e (right) for 2007. Only retrievals for QC = 0 and 1 are used. Shown are the global 
results as well as the day vs. night and the land vs. ocean distributions. 
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Fig. 13. Zonal averages of r (left column), 7c ice (center column), and D e (right column) for all of 2007, and also for DJF, MAM, JJA, and 
SON. The global results are in black (repeated from Fig. 12), while the seasonal day vs. night and land vs. ocean distributions are highlighted 
in color. 


further when placed in proximity to midlatitude cyclone cen- 
ters. 

4.1 Diurnal variations of ice clouds in the Maritime 
Continent 

The diurnal cycle of clouds, humidity, and precipitation has 
been quantified in recent years with global satellite data sets 


(e.g., Chen and Houze, 1997; Rossow and Schiffer, 1999; 
Yang and Slingo, 2001; Nesbitt and Zipser, 2003; Tian et al., 
2004, 2006). Although coarsely gridded traditional general 
circulation models (GCMs) continue to struggle in captur- 
ing the behavior of the diurnal cycle (e.g., Yang and Slingo, 
2001; Tian et al., 2004; Dai, 2006), the multiscale mod- 
eling framework (MMF or superparameterization) GCMs 
(e.g., Zhang et al., 2008; Pritchard and Somerville, 2009) or 
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global cloud-resolving models (e.g., Sato et al., 2009) have 
proven to be more representative of its amplitude, phase, and 
other complexities. The amplitude and phase of the diur- 
nal cycle are strongly dependent on the region (e.g., land- 
sea contrast, topography) and physical parameter of interest, 
including cloud-related quantities like precipitation, cloud 
fraction, height, and r (Cairns, 1995; Sato et al., 2009). 
A climate change signal in either the amplitude or phase 
of the diurnal cycle can have profound impacts on climate 
trends through the modulation of the daily timing of maxi- 
mum (or minimum) cloud reflection, and the absorption and 
re-emission of infrared radiation by high ice clouds (e.g., 
Cairns, 1995). Convective precipitation (Nesbitt and Zipser, 
2003; Dai, 2006), cloud frequency/amount (Chen and Houze, 
1997; Rossow and Schiffer, 1999; Tian et al., 2004), and out- 
going longwave radiation (e.g., Taylor, 2012) are perhaps the 
best-observed cloud-related quantities over the diurnal cycle. 

It has been well documented that there is a clear land- 
sea contrast for the diurnal cycle of high clouds (cloud tops 
above 440hPa) (e.g., Yang and Slingo, 2001; Tian et al., 
2004). High-cloud amount over tropical land is observed to 
have a distinct minimum during midday, with a maximum in 
the evening and at night (Cairns, 1995; Rossow and Schif- 
fer, 1999; Tian et al., 2004). In contrast, high-cloud amount 
over tropical oceans is observed to have a minimum during 
early morning and a maximum in the afternoon (e.g., Tian et 
al., 2004). Stubenrauch et al. (2006) shows that a maximum 
in tropical high-cloud amount is obtained over land during 
evening and night, while thin cirrus maximizes early in the 
afternoon, with an overall minimum in high-cloud amount 
around solar noon. Using CALIOP observations, Sassen et 
al. (2009) show that diurnal differences in thin versus thick 
ice cloud frequency strongly depend on the proximity to land 
and ocean, and can be of opposite sign depending on the 
range of r. 

Previous studies (e.g., Neale and Slingo, 2003; Tian et al., 
2006; Qian, 2008) have shown that the diurnal cycle is very 
strong in the Maritime Continent, and it may play a funda- 
mental role in the global climate. The AIRS instrument sam- 
ples two times during the diurnal cycle at 01:30 and 13:30 LT. 
There are geographical locations and geophysical parameters 
for which AIRS provides no diurnal variability information 
because maxima or minima occur in between the AIRS local 
sampling times. Fortunately, over the Maritime Continent, 
with its strong diurnal maxima and minima driven by heat- 
ing differences imposed by land-ocean contrasts, the local 
crossing time of AIRS lends itself well to sampling impor- 
tant aspects of the convective diurnal cycle. Initial results of 
the diurnal variability of AIRS ice cloud frequency, r, and 
D e are shown in Fig. 14. The diurnal cycle of ice cloud fre- 
quency is not especially large: generally less than 10-20%. 
In the case of r, the diurnal differences exceed a factor of 2 
and are highest over and adjacent to the islands. The min- 
imum (maximum) at 13:30 LT (01:30LT) over the islands 
(adjacent oceans) are consistent with the spatial and temporal 


Mean Value ASC-DES (1330-0130 LT) 



Fig. 14. AIRS diurnal variations of ice cloud properties over the 
Maritime Continent during 2007 for the mean value (left column) 
and the ascending-descending differences (13:30-01:30 LT, right 
column). Shown are ice cloud frequency (top row), r (middle row), 
and D e (bottom row). 


variations of TRMM convective features (Nesbitt and Zipser, 
2003). These patterns are also simulated by the 7 km NICAM 
(Sato et al., 2009), which shows evidence for a 3 h lag of the 
maximum in high-cloud amount behind the maximum in pre- 
cipitation (e.g., Tian et al., 2004). 

The diurnal differences for Z) e are less pronounced than 
for r, with some suggestion of structure near the individ- 
ual islands. However, D e appears to be spatially out of phase 
with r in many locations on a daily basis (not shown). The 
diurnal variability is also either more or less emphasized on 
either side of the traverse range of Papua New Guinea de- 
pending on the time period investigated (not shown). Protat 
et al. (2011) show that, using ground-based retrievals at the 
Darwin Atmospheric Radiation Measurement (ARM) pro- 
gram observing site, the magnitude of r e is dependent on 
the large-scale flow regime and the type of ice cloud present 
(deep convective, anvil cirrus, and thin cirrus). 

This initial snapshot of the diurnal cycle in the Maritime 
Continent confirms that its amplitude is dependent on cloud 
parameter and geographical location. The day-night differ- 
ences in AIRS demonstrate skill and are offering new in- 
sights on the diurnal variability of ice cloud properties, as 
well as thermodynamic phase in many other geophysical 
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regimes. Furthermore, FOV-scale matches of temperature, 
water vapor, and cloud properties can now be composited 
upwards from the FOV scale, preserving rich and detailed 
spatial and temporal variability. 

4.2 Midlatitude cyclone composites 

While there is ample observational evidence of a climate- 
change-induced poleward shift in the storm tracks in both 
hemispheres (Bengtsson et al., 2006; Johanson and Fu, 2009; 
Bender et al., 2012), change in the frequency and inten- 
sity of midlatitude cyclones within each of the storm tracks 
is much less certain (Schneider et al., 2010). One particu- 
larly successful and rigorous approach to evaluate present- 
day climate model simulations of midlatitude cyclones is 
to composite their dynamic and thermodynamic features in 
a common coordinate system relative to the cyclone center 
(e.g., Lau and Crane, 1995; Naud et al., 2006; Field and 
Wood, 2007). Field et al. (2008) used a number of CAM3 
physics perturbations to test this approach, including a new 
microphysics scheme, and their simulations were compared 
to cloud amount, rain rate, and near-surface winds derived 
from MODIS, AMSR-E, and QuikSCAT, respectively. Field 
et al. (2008) found that all model perturbations produced 
too much thick high cloud within cyclones, although impor- 
tant differences were found between the perturbations. The 
thick high-cloud bias was significantly reduced in the CAM5 
(Kay et al., 2012). CAM3 also showed strong relationships of 
high-cloud fraction with cyclone intensity and column water 
vapor amount, while satellite observations suggest a strong 
relationship only between high-cloud fraction and cyclone 
intensity (Field et al., 2008). Given the availability of tem- 
perature and water vapor mixing ratio profiles from AIRS in 
conjunction with the new cloud parameters reported in this 
work, this topic warrants continued investigation. 

Composites from the UK Met Office atmosphere model 
with ISCCP cloud type occurrence frequency and their short- 
wave radiative effects (Bodas-Salcedo et al., 2012) show 
that the SH surface shortwave bias in CMIP3 described by 
Trenberth and Fasullo (2010) is due largely to a dearth of 
low clouds in the cold sector of cyclones. This bias was 
partly mitigated by improving the representation of clouds in 
shear-dominated boundary layers. Flowever, Bodas-Salcedo 
et al. (2012) speculate that large biases may remain in an- 
ticyclones that also contain significant low cloudiness. Fur- 
thermore, biases remain in the cold air sector from poor sim- 
ulations of midlevel clouds; in particular the UKMO model 
places midlevel clouds too close to the cyclone center, and 
produces too few of them. 

The new collection of AIRS cloud products provides ad- 
ditional constraints for model evaluation. Using January 
2007 AIRS cloud property retrievals, we construct cyclone- 
centered composites separately for SH and NH oceanic cy- 
clones using the method of Naud et al. (2012). Despite the 
sample size limitations, AIRS composites of cloud cover are 


seen to be consistent with previous studies (e.g., Field and 
Wood, 2007) and the new products provide additional in- 
formation on clouds in cyclones. Here we discuss differ- 
ences and similarities between NH winter and SH summer 
cyclones. For convenience, SH cyclone composites are re- 
versed along the north-south axis, so that the poleward side 
is at the top of each figure. This allows direct comparison 
between the two hemispheres. 

The total ECF, upper layer 7c, and lower layer Tq fields for 
the NH and SH composites are shown in Fig. 15. The highest 
ECF occurs along the warm front in both hemispheres, while 
elevated values extend poleward and along the cold front in 
the SH (Field and Wood, 2007). The upper-level 7fc is cold- 
est to the north of the cyclone center in the NH. In the SH, 
the coldest upper-level Tq is found in the warm front region, 
and the warmest Tq is equatorward, at the back of the cold 
front, where open-cell convective tops generally occur at low 
altitudes. For lower-level 7c, a similar pattern to upper-level 
7c is seen in the NH (aside from higher values of 7c); how- 
ever, the coldest lower-level 7c is rotated poleward in the SH 
compared to the upper-level 7) . This suggests fundamental 
differences in multilayer cloud structure between SH and NH 
cyclones or between summer and winter cyclones. These fun- 
damental differences deserve a thorough study with multiple 
years of data. Recall that these particular cloud fields include 
all cloud types and structures. Below, the composites will be 
sorted by cloud thermodynamic phase. 

The cloud occurrence frequency of all cloud types and the 
three phase categories (ice, liquid, and unknown) are shown 
in Fig. 16. The cloud frequency (Fig. 16) and ECF (Fig. 15) 
fields have larger differences poleward of the cyclone center 
in the NH compared to the SH. This strongly suggests that 
clouds in the NH cyclone composites are optically thin pole- 
ward of the storm center. Ice clouds are most common pole- 
ward and eastward of the cyclone center in the NH (Field and 
Wood, 2007), and this pattern is rotated slightly more east- 
ward in the SH, with an extension equatorward into the warm 
sector not seen in the NH. However, the peak frequency is 
higher in the NH, as with total cloud frequency. There is a 
very high occurrence of liquid cloud poleward and west of 
the cyclone center in the SH, while an opposite pattern ap- 
pears in the NH with much less liquid frequency overall. The 
unknown category is most frequent in the cold sector behind 
and along the cold front and has a higher magnitude in the 
NH. This is consistent with the presence of open-cellular cu- 
mulus and the difficulty of assigning cloud phase because the 
cloud tops are frequently located in the mixed-phase temper- 
ature range (e.g., Nasiri and Kahn, 2008; Klein et al., 2009). 

The ice cloud properties are shown in Fig. 17, with a com- 
posite image of ECF restricted to ice clouds. The SH and NH 
patterns and magnitudes of ice cloud ECF are more similar 
than that for the total ECF shown in Fig. 15. There are some 
important differences between the SH and NH cyclones high- 
lighted by Fig. 17. The warm (cold) front is more prominent 
in the NH (SH). In both hemispheres, the highest values of 
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Fig. 15. Midlatitude cyclone composites for January 2007 for the total ECF, upper-level Tq, and lower-level Tq in the (a)-(c) NH and (d)-(f) 
SH. 
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Fig. 16. Midlatitude cyclone composites of cloud frequency and ice, liquid, and unknown phase clouds for (a)-(d) the NH and (e)-(h) SH 
during January 2007. 


ice cloud ECF are closer to the cyclone center when com- 
pared to the total ECF in Fig. 15. The r patterns also track 
ice cloud ECF patterns. There are higher values of r in the 
NH, and again a greater emphasis on the warm (cold) front 
in the NH (SH). In the case of D e , higher values occur equa- 
torward of the cyclone center in both hemispheres. However, 
higher values are found in the SH, and the hemispheric dif- 
ferences are larger poleward of the cyclone center. 7c, ice is 
much colder than the upper-level Tq shown in Fig. 15 and 
is consistent with it being an ice cloud subset of all cloud 
types. Ice clouds are slightly colder in the SH compared to 
the NH, which is surprising considering that these data are 


from the austral summer and boreal winter. This may indi- 
cate that large-scale cloud ice conditions are always colder 
in the SH, as found when comparing winter cyclones in both 
hemispheres in Naud et al. (2012). The coldest 7c, ice are 
found poleward and eastward of the cyclone center and the 
warmest are found in the cold sector in the NH (consistent 
with Fig. 15), although the relative frequency of ice in the 
cold sector is only between 10 and 30 % (Fig. 16). In the SH, 
this contrast is slightly rotated so most of the eastern side ice 
clouds are colder than their western side counterparts. Again, 
this could be a seasonal feature, which we intend to explore 
further as a large set of AIRS data becomes available. 
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Fig. 17. Same as Fig. 16, except for ice cloud ECF, r, D t , and 7c, ICE- 


5 Discussion and summary 

We describe the retrieval methodology and initial results of 
the Atmospheric Infrared Sounder (AIRS) version 6 (v6) and 
Advanced Microwave Sounding Unit (AMSU) (Chahine et 
al., 2006) instrument suite cloud products. The cloud top 
properties (temperature/pressure/height and effective cloud 
fraction) are obtained for up to two layers in the AIRS 
level 2 (L2) standard product, and are now reported at the 
AIRS field-of-view (FOV) resolution. Significant improve- 
ments in cloud height assignment over version 5 (v5) are 
shown with FOV-scale comparisons to cloud vertical struc- 
ture observed by the CloudSat 94 GFlz radar (Stephens et 
al., 2008) and the Cloud-Aerosol LIdar with Orthogonal Po- 
larization (CALIOP; Winker et al., 2010). These improve- 
ments are obtained for most observing conditions including 
land, ocean, day, and night. More realistic small-scale cloud 
structures are now observed, partly due to improvements in 
boundary layer cloud characterization. 

Several new cloud products are now reported in the AIRS 
L2 support product files. The first is the cloud thermody- 
namic phase that identifies clouds as “ice” or “liquid”, with 
an additional category “unknown” for confidently detected 
clouds that do not contain an easily identifiable ice or liq- 
uid signature. Jin and Nasiri (2013) showed that AIRS very 
consistently detects ice, but not liquid, when compared to 
CALIOP. Ice cloud frequencies in excess of 80-90% are 
found over the Maritime Continent, tropical South Amer- 
ica and central Africa. Larger amounts of ice are found in 
the winter storm track compared to the summer storm track. 
Very high frequencies of liquid cloud occurrence are detected 
around Antarctica in the austral summer and the patterns and 
magnitudes are consistent with previous studies (e.g., Hu et 
al., 2010; Choi et al., 2010; Wood, 2012). Although the stra- 
tocumulus regions contain primarily liquid phase, the trade 


cumulus regions are dominated by unknown phase, and this 
is consistent with the weak IR radiance cloud-phase signal 
from these cloud types. 

Three other new cloud products include ice cloud effec- 
tive diameter ( D e ), ice cloud optical thickness (r), and ice 
cloud top temperature (7c, ice), and are derived using an op- 
timal estimation approach adapted from the Tropospheric 
Emission Spectrometer (TES; Bowman et al., 2006) cloud 
retrieval methodology restricted to AIRS FOVs that contain 
ice clouds. Both r and D e are retrieved simultaneously along 
with 7c, ice to obtain better fits and more frequent retrieval 
convergence than if r and D e are retrieved alone. Quality 
control (QC) parameters are described that streamline the use 
of these ice cloud properties and are based on the quality of 
radiance fits between simulations and observations and the 
magnitudes of the averaging kernels. Distributions of r and 
D e for 2007 show that r is highest in the deep tropics and 
oceanic midlatitude storm tracks, and lowest in the subtrop- 
ics, the Arctic, and over Antarctica. The region of high r in 
the deep tropics is much more confined in latitude than the 
high occurrence frequency of ice cloud, which extends to the 
edge of the subtropics; this pattern captures the narrow region 
of deep convective clouds and adjacent thin cirrus. The distri- 
butions of D e peak between 30 and 40° N/S, with minima in 
the broad region of the tropics and high latitudes (Curry and 
Ebert, 1992). A small increase in D e associated with deep 
tropical convective events is shown to be consistent with both 
MODIS (Yuan and Li, 2010) and surface-based radar and li- 
dar retrievals (Protat et al., 201 1). 

The zonal mean differences of the diurnal variability of 
ice cloud properties are not as revealing as regional-scale 
differences. The regional context is important for grasping 
the rich complexity of the diurnal variability that is actu- 
ally present in AIRS data. The diurnal variability of r in the 
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Maritime Continent is shown to be significantly larger than 
diurnal variability for the total cloud fraction, ice cloud fre- 
quency, and D e , and is linked to the spatial pattern of the 
convective lifecycle that is anchored to the island archipelago 
morphology. The distributions and magnitudes are in sig- 
nificant agreement with seminal observations made by the 
Tropical Rainfall Measuring Mission (TRMM; Nesbitt and 
Zipser, 2003) and the International Satellite Cloud Climatol- 
ogy Project (1SCCP; Rossow et al., 2005) data. The infrared- 
based retrievals of AIRS are providing unique decadal-scale 
and global snapshots of clouds within the diurnal and annual 
cycles at all latitudes. 

The interhemispheric differences in the total cloud amount 
are rather small. However, the cloud-phase frequencies and 
ice cloud properties show substantially larger differences. To 
put the problem into an atmospheric process-based frame- 
work, we composited cloud properties using the approach of 
Naud et al. (2006) to quantify structural differences between 
winter NH and summer SH midlatitude cyclones. Although 
only one month (January 2007) is shown, which represents 
the boreal winter and austral summer, the observed struc- 
tural differences of cloudiness in proximity to the warn and 
cold fronts is consistent with previous works. This includes a 
very high frequency of liquid clouds poleward of summer SH 
cyclone centers, and a much larger occurrence of unknown 
clouds in open-cellular convection in the cold sector behind 
the winter NH cold fronts. The frequency of cirrus is high- 
est along the warm (cold) front in the NH (SH), consistent 
with Naud et al. (2012). The new ice cloud property retrievals 
show elevated D e along the warm and cold fronts, and much 
lower values poleward of the cyclone centers in both hemi- 
spheres. However, the overall values of the boreal winter D e 
are lower than in the austral summer. Further research will 
quantify the full range of seasonal and interannual variabil- 
ity. 


6 Future work 

Many improvements in the retrieval approach should be in- 
vestigated further. Multilayered clouds could be added to 
the scattering calculation. Since AIRS reports two layers in 
the standard retrieval, including the lower layer in a radia- 
tive transfer calculation is a possibility. However, the ap- 
proach for treating the lower (or additional) layer(s) raises 
a large set of complications. As phase is only assigned for 
the top of the upper cloud layer, assigning a phase for a 
lower layer is ambiguous (and untested), but could be ap- 
proximated by retrieved values of Tq. Additionally, it is not 
certain whether the lower layer will have a signal that is un- 
ambiguous enough to retrieve the optical and microphysical 
properties, whether it is liquid, ice, or possibly mixed phase. 
Another approach is to limit the most rigorous retrievals to 
the CloudSat/CALIPSO track and better constrain the verti- 


cal structure and phase, but it is unclear how much additional 
information the active sensors will provide. 

The ice cloud property retrieval is a post-processor that 
runs after the AIRS standard level 2 cloud-clearing algo- 
rithm. It is possible to include some (or all) elements of 
the cloud retrieval at earlier steps in the cloud-clearing pro- 
cess to improve the overall L2 geophysical retrieval results. 
For instance, a simple estimate of the cloud thermodynamic 
phase in the iterative cloud-clearing steps may improve the 
L2 full geophysical state retrieval. The record of collocated 
matchups of AIRS and MOD1S data at the FOV/pixel scale 
is now over 10 yr long. Advancements in the instrument cal- 
ibration, ground-truth comparisons through validation, and 
the collocation methodology have brought the two instru- 
ments closer to a seamless pan-spectral sensor that, in the- 
ory, could be used for joint retrievals. Better estimates of 
the a priori cloud structure from CloudSat, CALIPSO, and 
other satellite platforms will help improve the ice cloud prop- 
erty retrievals and may improve some of the thermodynamic 
phase estimates. The retrieval will benefit from improve- 
ments in other information, including prior covariances in the 
geophysical phase space that may depend strongly on cloud 
regime. 


Appendix A 


Acronym 

Meaning 

AIRS 

Atmospheric Infrared Sounder 

AMSR-E 

Advanced Microwave Scanning Ra- 
diometer - Earth Observing System 

AMSU 

Advanced Microwave Sounding Unit 

ARM 

Atmospheric Radiation Measurement 
program 

AVHRR 

Advanced Very High Resolution Ra- 
diometer 

CALIOP 

Cloud-Aerosol LIdar with Orthogonal 
Polarization 

CAM3/CAM5 

NCAR Community Atmospheric 
Model (CAM) version 3 or 5 

CMIP3/CMIP5 

Coupled Model Intercomparison 
Project Phase 3 or 5 

CrIS 

Cross-track Infrared Sounder 

D4S 

Delta-four-stream approximation for 
radiative transfer 

ECF 

Effective cloud fraction 

ECHAM5 

Max Planck Institute for Meteorology 
atmospheric general circulation model 

EUMETSAT 

European Organisation for the Ex- 
ploitation of Meteorological Satellites 

FOR 

Field of regard (3x 3 AIRS FOVs 
coregistered to AMSU FOV) 

FOV 

Field of view (AIRS footprint) 


www.atmos-chem-phys.net/14/399/2014/ 


Atmos. Chem. Phys., 14, 399-426, 2014 


420 


B. H. Kahn et al.: The Atmospheric Infrared Sounder version 6 cloud products 


HIRS 

High-resolution Infrared Radiation 
Sounder 

IASI 

Infrared Atmospheric Sounding Inter- 
ferometer 

1PCC AR4 

Intergovernmental Panel on Climate 
Change Fourth Assessment Report 

ISCCP 

International Satellite Cloud Climatol- 
ogy Project 

JPSS 

Joint Polar Satellite System 

L2 

Level 2 satellite data (retrieved prod- 
ucts at FOV scale) 

L3 

Level 3 satellite data (gridded FOV 
scale data for various time/space aver- 
aging) 

LBLRTM 

Line-By-Line Radiative Transfer 
Model 

MMF 

Multiscale Modeling Framework 

MODIS 

Moderate Resolution Imaging Spectro- 
radiometer 

NASA 

National Aeronautics and Space Ad- 
ministration 

NICAM 

Nonhydrostatic ICosahedral Atmo- 
spheric Model 

NOAA 

National Oceanic and Atmospheric Ad- 
ministration 

QC 

Quality control flag 

SARTA 

Stand-Alone AIRS Radiative Transfer 
Algorithm 

Suomi NPP 

Suomi National Polar-orbiting Partner- 
ship 

TES 

Tropospheric Emission Spectrometer 

TRMM 

Tropical Rainfall Measuring Mission 

TTL 

Tropical tropopause layer 

UKMO 

United Kingdom Meteorological Of- 


fice 


Appendix B 

Below is a simplified overview of cloud-clearing steps and 
effective cloud fraction/cloud top pressure retrievals. 

The AIRS retrieval is a cloud-clearing retrieval in which 
all parameters except clouds are retrieved from the “clear” 
radiances produced by the cloud-clearing process (Susskind 
et ah, 2003). Cloud-clearing selects the linear combination 
of the nine sets of observed radiances that best matches ra- 
diances calculated from the current estimate of the physical 
state, but without clouds, thus representing radiances from 
the clear portion of the field of regard (FOR). Once deter- 
mined, this same linear combination is applied to all channels 
except those above the clouds. 

The clouds are fit to the signal that is the observed radi- 
ances minus the calculated radiances from the current esti- 
mate of the clear state (Susskind et ah, 2006). This is done in 
an iterative manner with the following steps: 


2. Perform a microwave-only temperature and water va- 
por retrieval to adjust the first guess to best match ob- 
served microwave radiances. 

3. First cloud solution on AMSU FOR, starting from an 
arbitrary cloud state and minimizing (two cloud top 
pressures and two effective cloud fractions per AMSU 
FOR). 

4. First cloud-clearing solution. 

5. Minor adjustments to geophysical state. 

6. Second cloud solution on AMSU FOR, starting from 
step (3) cloud state (two cloud top pressures and two 
effective cloud fractions per AMSU FOR). 

7. Second cloud-clearing solution. 

8. Physical retrieval updating reflectivity, surface temper- 
ature and shortwave emissivity, the temperature and 
water vapor profiles, longwave emissivity, and the 
ozone profile, in that order. 

9. Third cloud solution, starting from step (6) cloud state. 
This one fixes the pressure on the AMSU scale but has 
separate cloud fractions on the AIRS scale (2 cloud top 
pressures and 18 effective cloud fractions per AMSU 
FOR). 

10. Third cloud-clearing solution. 

11. Physical retrieval updating surface temperature and 
shortwave emissivity, the temperature profile, long- 
wave emissivity, carbon monoxide profile, and 
methane profile, in that order. 

12. Final cloud solution, starting from the step (9) cloud 
state, solving for both pressure and cloud fraction on 
the AIRS FOV (18 cloud top pressures and 1 8 effective 
cloud fractions per AMSU FOR). 

13. Final cloud-clearing and OLR calculation for output. 

14. Cloud thermodynamic phase discussed in Sect. 3.1 
and ice cloud optical thickness and effective diameter 
cloud retrieval discussed in Sect. 3.2. 


1 . Generate first guess of clear state from MIT stochas- 
tic cloud-clearing and neural network solution, supple- 
mented by climatologies for trace gases, upper atmo- 
sphere temperature and water, and surface emissivity. 
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Appendix C 


Table Cl. A list of the AIRS channels used in the retrieval (out of a total of 2378). The noise-equivalent delta temperature (NEdT) is 
reported in Kelvin for a scene T\, of 250 K. The center of the channel spectral response function (SRF) and the full width at half max- 
imum (FWHM) of the SRF are both reported in wave numbers (cm -1 ). These data are taken from the AIRS channel properties file 
L2.chan_prop.2012.07.01.v9.5.3.anc, which is available at http://disc.sci.gsfc.nasa.gov/AIRS/documentation/v6_docs. 


Ret 

Chan # 

AIRS 
Chan # 

SRF Center 
Freq 

NeDT 

SRF 

FWHM 

Ret 

Chan # 

AIRS 
Chan # 

SRF Center 
Freq 

NeDT 

SRF 

FWHM 

1 

150 

692.76 

0.360 

0.480 

31 

730 

890.41 

0.142 

0.623 

2 

186 

702.74 

0.289 

0.495 

32 

757 

899.62 

0.128 

0.635 

3 

204 

707.85 

0.261 

0.503 

33 

764 

902.04 

0.193 

0.638 

4 

221 

712.74 

0.256 

0.510 

34 

799 

921.64 

0.161 

0.667 

5 

239 

717.99 

0.237 

0.518 

35 

807 

924.56 

0.183 

0.671 

6 

244 

719.47 

0.259 

0.520 

36 

812 

926.39 

0.133 

0.674 

7 

256 

723.03 

0.240 

0.526 

37 

819 

928.96 

0.101 

0.678 

8 

295 

734.15 

0.364 

0.558 

38 

824 

930.81 

0.149 

0.681 

9 

338 

747.60 

0.359 

0.580 

39 

833 

934.16 

0.145 

0.686 

10 

369 

757.61 

0.224 

0.595 

40 

839 

936.41 

0.188 

0.689 

11 

472 

800.01 

0.418 

0.667 

41 

844 

938.29 

0.148 

0.692 

12 

478 

802.19 

0.269 

0.671 

42 

848 

939.79 

0.137 

0.695 

13 

482 

803.65 

0.302 

0.674 

43 

853 

941.69 

0.123 

0.698 

14 

503 

811.41 

0.283 

0.687 

44 

860 

944.35 

0.117 

0.702 

15 

510 

814.03 

0.415 

0.692 

45 

869 

947.80 

0.141 

0.707 

16 

516 

816.29 

0.514 

0.696 

46 

873 

949.34 

0.147 

0.709 

17 

521 

818.18 

0.253 

0.699 

47 

882 

952.83 

0.155 

0.715 

18 

530 

821.60 

0.384 

0.705 

48 

897 

958.69 

0.137 

0.724 

19 

541 

825.82 

0.315 

0.713 

49 

908 

963.04 

0.120 

0.730 

20 

549 

828.91 

0.431 

0.718 

50 

914 

965.43 

0.121 

0.734 

21 

553 

830.47 

0.304 

0.721 

51 

1181 

1093.94 

0.146 

0.968 

22 

565 

835.17 

0.397 

0.729 

52 

1195 

1101.11 

0.612 

0.981 

23 

577 

839.92 

0.313 

0.737 

53 

1200 

1103.69 

0.164 

0.985 

24 

585 

843.11 

0.424 

0.743 

54 

1211 

1109.40 

0.233 

0.996 

25 

653 

865.21 

0.217 

0.590 

55 

1218 

1113.07 

0.167 

1.003 

26 

667 

869.68 

0.209 

0.596 

56 

1223 

1115.70 

0.249 

1.007 

27 

687 

876.15 

0.280 

0.604 

57 

1235 

1122.06 

0.338 

1.019 

28 

695 

878.77 

0.179 

0.608 

58 

1247 

1128.50 

0.210 

1.031 

29 

711 

884.05 

0.208 

0.614 

59 

1257 

1133.91 

0.249 

1.041 

30 

722 

887.72 

0.141 

0.619 
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